Palatini f{R) Black Holes in Nonlinear Electrodynamics 



O 

> 

o 
O 



> 

o 

00 

o 
d 



X 



Gonzalo J. Olmo^ and D. Rubiera-Garcia^ 
^ Departamento de Fisica Teorica and IFIC, Centra Mixta Universidad de 
Valencia - CSIC. Umversidad de Valencia, Burjassot-46100, Valencia, Spain and 
^ Departamenta de Fisica, Universidad de Oviedo, 
Avenida Calva Satelo 18, 33007, Oviedo, Asturias, Spain 
(Dated: October 6, 2011) 

The electrically charged Born-Infeld black holes in the Palatini formalism for f{R) theories are 
analyzed. Specifically we study those supported by a theory f{R) = R ± R? /Rp, where Rp is 
Planck's curvature. These black holes only differ from their General Relativity counterparts very 
close to the center, but may give rise to different geometrical structures in terms of inner horizons. 
The nature and strength of the central singularities are also significantly affected. In particular, for 
the model f{R) = R — R^ /Rp the singularity is shifted to a finite radius, r+, and the Kretschmann 
scalar diverges only as l/(r — r-+)^ . 
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I. INTRODUCTION 

Black holes are one of the most intriguing objects of 
Nature. Their event horizons dramatically modify the 
causal structure of space-time generating regions from 
where nothing can scape classically. From a semiclas- 
sical perspective, however, event horizons are responsi- 
ble for the generation of quantum radiation that slowly 
forces their evaporation in a manner already described 
by Hawking almost 40 years ago jH-Q. In terms of the 
Planck length, Ip = ^ hG/c^ ~ 10~'^^m, the event hori- 
zon of macroscopic black holes is located very far away 
from the center, where a singularity is expected classi- 
cally. This singularity, however, should be an artifact of 
the classical theory, and a full quantum theory of grav- 
ity should provide suitable mechanisms to avoid its for- 
mation. Within such a theory, the weak, perturbative 
quantum phenomenon of Hawking radiation should be 
consistent with the strongly non-perturbative quantum 
processes expected to happen at the black hole center to 
remove the classically-expected singularity. Therefore, 
in the vicinity of the singularity, the internal structure 
of black holes must be different from that predicted by 
classical general relativity (GR). 

In this paper we address the issue of the internal struc- 
ture of black holes from a purely classical perspective 
by considering an f{R) gravitational action containing a 
correction term that generates non-perturbative effects 
at the Planck scale. This is achieved by formulating the 
theory a la Palatini i.e., by considering metric and 
connection as independent fields. As a result, the dy- 
namical degrees of freedom of the theory are the same 
as in GR, which implies that there is no multiplicity of 
solutions. In fact, the number of solutions is the same as 
in GR, and their structure is significantly deformed near 
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(and only near) the characteristic scale set by the cor- 
recting terms of the action. This means that the internal 
structure of these black holes is different from that of GR 
only near the Planck scale, while for lengths much larger 
than Ip the departures from GR are totally negligible. 
We hope that the analysis of this type of solutions allows 
to provide new insights on how quantum gravitational 
effects could affect the interior of black holes. 

The modified dynamics generated by the Palatini ver- 
sion of f{R) theories is intrinsically different from that 
of their metric counterpart, where the connection is as- 
sumed to be defined in terms of the Christoffel symbols 
of the metric. In the usual (metric) formulation, the 
modified dynamics is due to the existence of an effective 
dynamical scalar degree of freedom that can be defined 
as = df/dR. In the Palatini approach, though a scalar- 
tensor representation is also possible, the resulting scalar 
is non-dynamical and, therefore, there are no new degrees 
of freedom [5]. The equation for the independent con- 
nection can be solved by introducing an auxiliary metric 
that is conformally related with the physical metric g^i,. 
As we will show later, the conformal factor turns out to 
be a function of the trace, T = g^^^T^^, of the energy- 
momentum tensor T^i, of the matter fields. As a result, 
the field equations for the metric contain a number of 
new terms that explicitly depend on T and its deriva- 
tives. When r = 0, the field equations boil down to 
those of GR (with possibly an effective cosmological con- 
stant, depending on the choice of f{R) Lagrangian), but 
when r 7^ 0, the connection-induced matter terms yield 
modified dynamics 3]. According to this. Palatini /(i?) 
theories have the same (classical) vacuum solutions as 
GR with a cosmological constant. However, if semiclas- 
sical effects are taken into account (Hawking radiation, 
trace anomaly, . . . ) then black hole solutions should be 
different for different /(i?) Lagrangians because then the 
trace T would no longer be zero. 

Besides semiclassical scenarios, which are difficult to 
treat consistently, non-vacuum spacetimes that are of 
interest to explore Planck-scale effects of the Palatini 
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dynamics are those of electrically charged black holes. 
Though the stress-energy tensor of Maxwell electrody- 
namics is traceless and, therefore, cannot be used to 
probe the Palatini dynamics, it seems reasonable to ex- 
pect that besides the gravitational field, other fields could 
develop non-linear corrections at high energy scales. 
Such is the case of nonlinear electrodynamics (NED), 
which are expected to arise as modifications of Maxwell 
theory in several scenarios, such as effective Lagrangians 
in Quantum Electrodynamics @. Moreover, in order to 
break the tracelessness of the energy-momentum tensor 
and obtain electrically charged non-standard black holes 
NEDs is the natural choice. In this sense a physically well 
motivated candidate to be studied is Born-Infeld theory. 
The Born-Infcld (BI) Lagrangian was introduced in the 
thirties 7] in order to remove the divergence of the elec- 
tron's self-energy in Classical Electrodynamics. More- 
over, BI theory is singled out among the class of NEDs by 
its special properties concerning wave propagation such 
as the absence of shock waves and birefringence phenom- 
ena [8|], and in addition it enjoys an electric-magnetic du- 
ality [9|. In modern times, BI theory has received special 
attention as it arises in the low-energy regime of string 
and D-Brane physics [loj . 

Electrically charged BI black hole solutions within the 
context of asymptotically flat Einstein gravity has been 
thoroughly studied and characterized for many years (see 
e-g- uM)- Other examples of gravitating non-linear elec- 
trodynamics have been considered [12I], and moreover a 
general analysis on lar ge f amilies of gravitating NEDs has 
been also carried out [IJ] . Extensions of BI black holes 
to asyrnptotically (Anti-) de Sitter spaces can be found 
in Ref.[ij]. Let us finally mention that NEDs have also 
been used as attempts to classically remove the black hole 
curvature singularity [l5| . 

Beyond GR, black hole solutions within the context 
of higher order gravity theories (e.g. Lovelock theory 
[l6j ) have attracted much attention in the last years, as 
these theories are suggested to arise as corrections to 
the Einstein-Hilbert action by some low-energy results 
of string theory [ITj . In this context some authors have 
looked for black holes in Gauss-Bonnet theory coupled 
to BI electrodynamics [l8|, as well as to related theories 
[Ti^] , finding large differences in terms of horizons and sin- 
gularities as compared to their GR counterparts as well 
as striking new features such as the existence of branch 
singularities [20| . 

Literature on black hole solutions in the Palatini ap- 
proach is rather scarce and has focused essentially on 
vacuum solutions of f{R) theories. Most research at this 
regard deals with thermodynamics issues (see e^. [2lj 
and references therein) and entropy definition |22| . The 
aim of this paper is to provide an analysis on the ge- 
ometrical structure of electrically charged black holes, 
through NEDs, and to compare it to the GR limit. In 
particular, we look for the deviance from GR of the struc- 
ture of the BI black holes near their center, where the 
effects introduced by the Planck-scale corrected model 



f{R) = R zt R^/Rp become relevant (throughout the 
paper Rp = Ij,^). These models have been studied pre- 
viously in the cosmological context, where it has been 
shown that they can avoid the big bang singularity in 
isotropic and homogeneous scenarios [23i - i25i] , and also 
in astrophysical [26j . laboratory [23|, and solar system 
settings [2^. We provide a thorough analysis and dis- 
cussion of the metric components, Kretschmann scalar, 
and inner horizons and characterize them according to 
the values of the relevant length scales involved. We find 
that, in general, the metric components can be solved in 
terms of power series expansions, which in some partic- 
ular cases can be explicitly written in terms of special 
functions. Though the central singularity still appears in 
these models, we find that for the case f{R) = R—R^/Rp 
it is shifted to a finite radius, which we denote r+, and 
its strength as measured by the Kretschmann scalar is 
significantly softened as compared to GR or the model 
f{R) — R + R^ /Rp, for which it occurs at r = 0. 

The paper is organized as follows. In section II we 
introduce the general setup and basic elements for both 
the modified gravity /(i?) and matter NED source, and 
solve the Einstein equations in Palatini formalism. In 
section III we introduce the BI model, that shall be 
used throughout the paper. Section IV is devoted to 
the study of black holes within f{R) — R + R^ / Rp grav- 
ity theory coupled to BI NED, through the analysis of 
the Kretschmann scalar and the existence and features 
of the inner horizons. Then we extend this analysis to 
the f{R) = R — R^/Rp theory in section V and we con- 
clude in section VI with a summary and some future 
perspectives. 



II. DEFINITIONS AND FIELD EQUATIONS 

The action for GR coupled to Maxwell field is 

= 2^/ '^'''^^'^ J d^^^TgF^pF'^P , (1) 
87rG/c3, R = g^'^R^.iV), R^,{r) = RP^p,, 



-^.r^^+r^.r^^-r^^^r;;. represents the 
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R"" PtJ-v 

components of the Riemann tensor, the field strength of 
the connection P^^, and F^^^ = df^A^ — d^A^ is the field 
strength tensor of the electromagnetic vector potential 
A^. This definition leads to 
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Introducing the gauge field invariants as A" = 
-^FapF'^f' ^E^ -B^ and Y = -\FapF*°'f^ = 2E-B, 
with F*'^" = ^eP'^^'^Fap being the dual of the field 
strength tensor and E and B the electric and magnetic 
fields, respectively, the extension of (H]) to the f{R) case 
with NED source is just 



-g^{X,Y) . (3) 
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where ip{X,Y) is a given function of the two field in- 
variants, and defines the particular NED model. Note 
that the condition ip{X,Y) = ip{X - Y) for the NED 
Lagrangian density must be satisfied in order to im- 
plement parity invariance. Variation of the action ([3]) 
with respect to the metric for a radial electric field 
(£;(r)f ^0,B = 0), which implies that F = 0, yields 
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where we have introduced the notation fn = df / dR and 
tpx = §x- Taking the trace of (jlj with the metric g^"^ , 
we find that the scalar R is algebraically related to the 
energy-momentum tensor through 



RJr - 2/ 



(5) 



theory becomes equivalent to GR-(-A up to a constant 
rescaling of units. 

The electromagnetic field satisfies the equation 
{ifxF'"' + ^yF*'"') = 0. For a diagonal metric with 
spherical symmetry and a purely radial electric field, for 
which only F*^ is not zero, this equation becomes 
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(10) 

where q is an integration constant identified as the elec- 
tric charge for a given model. From this expression for 
F*-^ and the fact that X = —gttgrr{F*^^Y , it follows that 
for any spherically symmetric metric 
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This algebraic equation generalizes the GR relation R = 
— K^T to non- linear Lagrangians and its solutions will be 
denoted hy R = R{T). 

The variation of ([3]) with respect to F^^^ must vanish 
independently of (|4]) and gives (we assume a torsionless 
connection for simplicity [J|) 



-gfRg'"") = , 
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With this equation once a <p{X, Y — 0) theory is pro- 
vided, the solution for X = X(r) can be found, in prin- 
ciple, algebraically. 

To continue with our analysis, it is convenient to ex- 
press the NED stress-energy tensor as 
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where fa = fii{R[T]) must be seen as a function of the 
matter by virtue of Q. This equation can be solved 
easily (see for instance [231) ^^i^ leads to 



(7) 



F^^ is the Fevi-Civita connec- 



where /i^ 

tion of a metric h^^ conformally related with g^i,, being 
the conformal factor the function = /fl(i?[T]). Using 
this expression for the connection, we can write Q as an 
Einstein-like equation for the physical metric g^^, 
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One should note that all the i?, /, and /_r terms on the 
right hand side of this system of equations are functions 
of the trace T of the matter energy-momentum tensor. 
According to this, if T = 0, those equations boil down to 
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where Ae// = {RofRg — /o)/2/i?,, is evaluated at the con- 
stant value Ro = R[T — 0] and plays the role of an 
effective cosmological constant. This observation high- 
lights the fact that the modified dynamics of Palatini 
f{R) theories can only arise if T 7^ 0. For T = 0, the 



Using matrix notation, with / representing the 2x2 iden- 
tity matrix and a 2 x 2 zero matrix, we find that 
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which can be used to write T,,'^ as 
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In order to guarantee the positive definiteness of the en- 
ergy density for any NED model the condition 

P^Tl = ^ (2^xE^ - V{X, r = 0)) > 0, (15) 

will be assumed in what follows. This positive energy 
condition guarantees the single branched and monotoni- 
cally decreasing character of the associated E{r, q) field, 
as a consequence of Ea. ([TT|) (see Ref.j^). 
It is easy to see that the trace of the object p4|) leads 

to 



(16) 



which vanishes for Maxwell electrodynamics f{X) = X. 
This observation is very important because, as mentioned 
above, the modified dynamics of Palatini f{R) theories 
can only arise if the trace T ^ 0. Maxwell's linear elec- 
trodynamics times a constant is the only one (for van- 
ishing Y) satisfying the T = condition. Any nonlinear 
generalization circumvents the tracelessness condition of 
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the energy-momentum tensor, thus providing a deviance 
within Palatini f{R) theories from their GR counter- 
parts, as we shall see at once. Note that when the Y 
invariant is considered, the above trace expression is ex- 
tended as 



(17) 



The set of families with vanishing trace in this case is 
formed by all conic surfaces in the {X, Y, ip) plane having 
the origin as a vertex; for example the set of planes of 
the form (p{X, Y) — aX + bY (a, b constants) and all 
of them lead to the same solutions as for the GR case. 
However, here we shall be only interested in the purely 
electric case. 

A glance at Egs. p^ and dTTT) indicates that for a given 
NED the dependence of T^" on the radial coordinate r 
can be completely known before solving explicitly the 
equations for the metric. With the knowledge of T^"^ one 
can proceed to solve for ij. To do it, one may con- 
sider directly ([8]), whose right hand side involves several 
derivatives of T, or may consider the equations for ft^^ 
given in ^ and then use the conformal transformation 
to obtain g^^^. For notational and technical simplicity, 
we choose the second option. We begin by writing (|4]) as 



which in matrix notation becomes 



(18) 







JR 

^" ' (19) 

We now need to find the differential equations for the 
metric /i^j^ and use the conformal relation between met- 
rics to obtain g^j^^. To simplify the field equations for h^^ 
we define the line element in the two frames as follows 



-gttdt^ + grrdr"^ 
1- (-A{f)e^^'^dt^ + 

jR \ 



r^dfl'^ = 
dr^ 



1 

Tr 



-ds 



^■2 



A{f) 



(20) 



Therefore, to solve for the metric g^i, we will solve first 
for hf^iy and then use the relation d^U]) to find f/^i,. Note 
in this sense that r and r are related by = /r. 

With this decomposition of and that choice of co- 
ordinate f , we can use the expressions for the components 
of Rf/{h) to obtain the following expressions 



— Rf^ — — ; — — 
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(21) 
(22) 



The first of these equations implies that ip{f) is a con- 
stant, which can be absorbed into a redefinition of the 



time coordinate and will thus be omitted. On the other 
hand, taking into account the relation = "t^Jr and 
defining A{f) = 1 — 2M (f) /f, we find that 



^Ir 



which in terms of the r variable becomes 



Mr 



providing the right expression in the GR hmit {fji — )■ 1) 
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The equation is far too complicated to be solved for 
a generic f{R) gravity and ip{X,Y = 0) NED model. 
Hence we are led to consider particular cases of theories 
for both the matter and gravity fields, as described in the 
following sections. 



III. BORN-INFELD NED 



In this case, the NED Lagrangian is defined as 



(26) 



For this theory and considering electrically charged solu- 
tions, Eg.pTj) becomes 



X 



1-^ r4 



which leads to 
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g2 + r4^2 



(27) 
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Since we will be dealing with expressions of the form 
it seems convenient to write that expression as 

follows 



1 



1 



X_ _ g2 ^ ^ 



where = With this notation we can write 



(29) 



^ = 2/32 1 - 



(30) 
(31) 
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where the limit to the usual electrodynamics corresponds 
to z > 1. 

For further reference we also need the expression for 
the BI Tj component 



- 1 



(32) 



which obviously satisfies the positive energy condition 
([T5|) for any /3. The energy associated to the ESS field is 
given by 



e((Z)=47r/ drr^T}{r,q)^ (33) 

r°o _3/2ol/2 

Indeed this quantity plays a relevant role in the charac- 
terization of the GR BI black holes [llj. Note that us- 
ing the scaling with q in Eg.lfM]). the quantity ^^^^jl^ = 

srls/i}'^ ~ 1-236 becomes a universal constant for a given 
parameter /3 in the BI model. 



more transparent if we consider the dimensionless vari- 
able M = M/{j^k'^I3^), which leads to 



(36) 



With this choice of variables, the only adjustable param- 



eter is A = , which appears in / and 



A. External horizon and the function G{z) 

Horizons arise when gu = 0. According to the de- 
composition of the metric that we have considered, this 
implies that gtt = —A{f)/fi^ = 0, which is equivalent to 
A{r) = (and also to Jr — > oo, which occurs at z = 0). 
According to our analysis, we have 



A{f) = 1 - 



2M(f) 



(37) 



where M(f) can be expressed as 



IV. STUDY OF f{R) =R + R'/Rp 

We now consider the quadratic model f{R) — R + 
R'^/Rp, for which R = -k^T, as follows from Eq.©. 
From the general expression (jl6l) we find that 
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, (34) 
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which can be expanded for large z as T//?^ 
O and for z ^ as T/jS^ « + 

If — h OlzY . This behavior at small z indicates that 

167r L J 

Jb — 1 — 2k^T/ Rp may vanish at some small z depend- 
ing on the sign and magnitude of the combination A = 
K^/3^/Rp. In particular, for f{R) = R + R^/Rp, Jr does 
not vanish anywhere, whereas for f{R) = R — R^ / Rp the 
function fp vanishes at a finite radius, which requires an 
independent analysis (see Section IV)) . 

In terms of the variable z, (l24l) can be written as 



R,z 



(35) 



where 7 = \/q/l3. Now we just need to integrate this 
function numerically in each case. Additionally, it is use- 
ful to note that the dimensions of the term (^f + j;^ 

.-2 fl2 



7V> 



are given by a factor k^/?"^, in other words, ^/ -|- ^V'J 

^'f3' (/ + 4^^) , where f ^ -f + f = T / 

and (p = f/ 0^ ■ We thus find that the analysis becomes 
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(39) 



where Mq is an integration constant that represents the 
Schwarzschild mass in the vacuum case. Since f = rfj^'^, 
the horizon condition becomes 



1-^M(.)^0 



(40) 



Using the definitions 7 = \fqf0, = 1/(k^/3^), r| = 
K^g^/(47r), and rs = 2A/o, the condition A{r) ~ be- 
comes 



l-|-2(47r) 




^G(z) 




/3 ^/-i/a 
(41) 

The exterior horizon occurs for z 3> 1, where jp k, 1 
and G[z) « 0. For astrophysical black holes (with low 
charge to mass ratio, Vq/rs <^ 1), it is easy to see that in 
this limit the left-hand side of (|^T|) is almost unity, which 
means that this horizon is located very near the point 

Since z = r/7, it follows 



z, « 1/ [(4^)1/4 (l^) 



that r/i = 7Z/j ~ rg, as expected. 

It should be noted that our way of writing (j4T|) high- 
lights the three scales involved in the problem, namely, 
the ratio charge-to-mass Vq/rs, the ratio NED-to-charge 
1/3 /rq, and the ratio NED-to-Planck given by A = 
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B. Solving for Af(z) 

A closed analytical expression for M{z) for arbitrary A 
is not possible in general. However, a glance at the series 
expansion of ([55)1 near z « (for finite A > 0), 

(TT-2X) (tt^ - 47rA + lOA^) „ , , , 
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327r3/2A3/2 



suggests that the choice A = 7r/2 may lead to some simpli- 
fications. Luckily this is indeed the case. For that choice 
of A we can find analytical solutions for that equation. 
For A = tt/2, Eq.dMl) becomes 



(3 + 2z4) (1 + (1 + z^) - iz^^l + z*) 

leV^TT (1 + Z*)^/-* (1 + 2z4)3/2 



(43) 

This expression can be compared with that of GR for the 
BI model 



1 + z4 - z^yrr? 



SvrVl + z4 



(44) 



which is simply —jp^ times Eq.(j32]), as expected from 
Eg.ipSj) and the units employed. In the BI case, inte- 
grating (|44)) we find that 



G{z)bi = — 
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iArcSinh 



(-l)iz 



-1 



(x/TT^-z2) 

247r 



127r 
Cgr 



(45) 



where the integration constant Cgr is necessary to get 
the right asymptotic behavior at z — > oo. Its value is 



Cgr — — 



2(-l)i/4EllipticK[-l] - (-l)3/4EllipticK[2] 



rl/2 



24r[3/4]2 



127r 

-0.0491809 



(46) 



Note that the physical interpretation of this constant is 
immediate, as it turns out to be — I/Stt times the univer- 
sal energy constant defined below Ea. (|33l) . 

In the /(i?) case with A = n/2, the solution is more 
complicated. In terms of a series expansion, for small z 
we find 
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P,{z) + P.(z)AFl [1, 1, -1, 1, -2z41 

48\/27r(l-f 3z4-f 2z8) 
(l + z4)i/4yTT^ [1 + 4^4(1 + ^4) 



- 12z2(2 + 3z'*)Vl + z4 
P2(z) = 24z2(l-F3z'' + 2z8) , 



(48) 
(49) 



where AF\ is the Appell hypergeometric function. In 
order to find a solution for z — >■ oo, it is convenient to ex- 
press p6p in terms of the variable x = 1/z. The solution 



for large z can then be expressed as 

_ Qi(z)-Q2(2)AF1 [\,=i,\ 
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48y27rz2 (1 -f 3z4 + 2z8) 
i(z) = (1 + z4) yrT2;J [z2 (1 + 4z4 (1 + z^)) 
4(4 



(50) 
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Q2{z) = 16V2z (l + 3z'^ + 2z^ 



(51) 
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These expressions for Gi^" and Gi^°° must be sup- 

2 2 

plemented with an integration constant. In the case of 
Gi^°°, the expansion at z — oo leads to Gi"^°° ~ 

2 2 

— -nr 1- ,nn^ B -I- O [i] , which recovers the usual GR 

expression at lowest order and, therefore, does not need 
the addition of any constant. The integration constant 
needed by G|r^° can be obtained by forcing the agree- 
ment between the G's in the overlapping region. Choos- 
ing any point where the two series are well defined, we 
find that G- = G^^'^(zo) - C^£^°(zo) « -0.02936. If 

2 2 2 

we take, for instance, zq — 1/a/2, we find 
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which fits with the value obtained above. Therefore, near 
the origin the solution is given by 



G? 



Gil -I- G! 



(54) 



The expansion of G near the origin leads to the follow- 
ing expression: 



GV° = Gn 



1 



3z4 



48V^7r 64V27r ' 
which suggests the definition of a new constant 
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in terms of which G^ 
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(55) 



(56) 
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This def- 



inition nicely fits with the expansion corresponding to 
general A, which is of the form 



Ga 



(7r-2A)z2 (^2 _ ^ 10A2) z4 
327r3/2VA ^ 



1287r3/2A3/2 



(57) 

It should be noted that the constant of integration C\ is 
negative from A = to A « 18.6444. 

As we will see at once, the A = 7r/2 solution is very 
useful to understand how the Planck length modifies the 
BI geometry of GR. Since we do not have complete an- 
alytical solutions for general A 7^ 7r/2, each case must 
be computed separately and, in this sense, the graphical 
comparison with that solution resulting from a numerical 
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G[z] 




0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 1. From bottom to top representation of the function 
G{z) for small and medium values oi X — j^,^, -I, tt, 27r, 47r 
(dashed curves). The solid green curve represents the special 
case A = -I and the solid blue one gives the GR (A — 0) case. 



G[z] 
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' , A=500;r 
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-- A=507r 
A=57r 
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-----2,6- 
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— j-.o 




A=;r/2 ^^^^S.:^—^ 










0.05 













FIG. 2. From bottom to top representation of the function 
G{z) for large values of A = 10 x f , 10^ x f , 10^ x f (dashed 
curves). The solid green curve represents the case X = ^ and 
the solid blue one the GR (A = 0) solution. 



integration will be very illustrative (see FigslHandE]). In 
this sense, it is important to note that the change in the 
function G\{z) observed in Fig [1] is smooth as A is con- 
tinuously increased from zero to positive values. This, 
however, requires an explanation, because from the ex- 
pansion one cannot recover ([^^ in the limit A — > 0. 
In (j42p we considered an expansion for z — >■ with finite 
A > 0. If we consider instead an expansion for A — >■ 
with finite z > 0, to first order in A we obtain 



i?a/375^"'^^* is given by 



Kret(r)^ + (59) 

where we are using the notation Ar = dA/dr. Given our 
decomposition of the metric, we have 



« Mf ^ -I- 



Ll(z) +i2(2)VT + ^ 



167r2z2 (H- z4) 



4x3/2 



A + 



(58) 



where Li(z) = z^-bz^-Az^^ and L2(z) = -2 + 3z''-|-4z^. 
This expression is in agreement with (j44p when A — >■ 
and is valid as long as the Planck scale corrections can be 
treated as small perturbations (recall that A = Ip/l'^) in 
the BI background of GR. The perturbative expansion in 
A, however, is not valid when we are near the singularity 
and the effects of the modified gravitational dynamics are 
important. Therefore, a faithful description of the geom- 
etry for small (but finite) A in the limit z — would 
require to sum over all the terms in the infinite pertur- 
bative series expansion ([55)1 or, equivalently, consider the 
non-perturbative expansion about z — >■ given in (1421) . 
This will be clearly seen in the following sections, where 
we study the properties of the singularity by computing 
the Kretschmann scalar. 



C. Kretschmann scalar 

Given a line element of the form ds"^ = —A{r)dt'^ 
B{r)dr'^ + r'^dft^, the Kretschmann scalar Kret 



Aiz) 



B{z) 



Jr 



(1 + (5i G'(z)) 



1/2-. 



1 - 



*2 2/r 



dz 



(60) 



(61) 



where the constants Si and 62 take the form 



61 = 2(47r) 




(62) 
(63) 



Since our solutions for the metric are series expansions, 
there is no chance to write the Kretschmann scalar in a 
simple form. We will content ourselves with series ex- 
pansions about the origin to check whether the singu- 
lar behavior of the metric improves or worsens in these 
models with respect to the GR solution. Indeed this is 
enough, as the f{R) theory chosen only deviates from the 
GR result precisely in that region. We will thus begin by 
computing the behavior of the BI model in the case of 
GR. 
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Bom-Infeld in GR 



In this case, an expansion around the point z — leads 



to 



Kret(z) 



+ 



(64) 



where Kret(z) = 7'*Kret(r), and the constant Cgr has 
been aheady defined in Ea. (l46l) . It is worth noting that 
for Maxwell electrodynamics the Kretschmann scalar is 
given by 



Kret(r) 



48Mo2 



48Afnr2 14rf 



+ 



(65) 



A[z] 




FIG. 3. Behaviour of the metric function A[z) for GR-BI 
black holes, according to the sign of 1 + CgrSi or, alterna- 
tively, of Mo — K'e{q). We are plotting the case &2 = 
Note that the other cases for Mq — K^e{q) may be obtained 
by moving the associated curve above or below the A — axis 
according to whether S2 > oi 62 < respectively. 



which diverges as ^ near the origin and indicates 
that the charge increases the intensity of the divergence. 
In the BI case, the maximum divergence is of order 
~ in general but if the combination of constants Si 
is constrained to take the value 1 + 5iCgr = (because 
Cgr — C\=Q < 0, see FigH]) then the divergence is soft- 
ened to the order of ~ 



-4. 



Kret(z) 



1 



1 



2 
3^ 



13 
15 



(66) 



Note that, as pointed out at the end of section pVBp . 
the series expansions (|64| and (|66l) are valid for finite 
values of Z^, while the limit Z^ — )■ would require to sum 
over the infinite terms of the series. This justifies why it 
is not obvious how to get (|55|) from taking the limit 
IjS ~> 0. Similar situations will arise in the A ^ cases 
discussed below. 

The function A{z) near the origin behaves as 



Aiz) 



1 + CguSi 
ZS2 



1 - 



Si 
SttSo 



z^Si z^Sx 



247r(52 8 7r(5, 



With the conditions 1 
this function becomes 



Cgr,'5i = and S2 



(67) 

1 



, ^2 z4 ^8 

Aiz) « \ 

^ ' 3 10 72 



(68) 



These special values of Si and S2, which are equivalent 
to constraining the ratios Vq/rs and rq/lp, receive a nice 
interpretation. Indeed, using the definitions introduced 
so far the first condition translates into 



Mo 



2 3/2 0I/2 



(69) 



so taking into account the value of Cgr given by pS)) 
and the expression of the BI field energy p3p this simply 



means Mq — K^e. On the other hand the second con- 
dition above may be written as n^q/S = 47r. These con- 
ditions represent limiting values of the BI and the black 
hole parameters splitting two different regions of black 
hole configurations (see Refs.flTI. [isjb The first condi- 
tion splits the metrics into those diverging to ±00 at the 
center, with equation (|69l) being the transition value be- 
tween both regimens, and for which the metric becomes 
finite at the center. In the latter case, the addition of the 
second condition makes the metric to vanish at the cen- 
ter, as shown in Eq. (|68ll . When Mq — K^e{q) < black 
holes with two horizons (inner and event), extreme black 
holes or naked singularities may be found (see Figl3]), 
while for Mq — K'^e{q) > a single event horizon appears. 

When 1 + CgrSi 7^ 0, the function B{z) near the origin 
behaves as 



Biz) « - 



S2Z 



62 {61 - 87r^2) 



With the conditions 1 
this function becomes 



1 + Cgr(5i 87r(l-|-CGR(5i) 
S2 {Si - 87r^2)^ 
64^2 (1 + CgkSi)^ ^ ■ ■ ■ ■ 

= and S2 — 



(70) 



SttCgr 



, , 3 9 



27z" 



507z6 



100 250 10000 



0[z]'. (71) 



Note that to obtain this last result one must impose the 
conditions on Si and S2 before expanding (|6ip about z ~ 
0. 



2. Bom-Infeld for A / 0, tv/2 



In this case the analytical expansion leads to 



Kret(z) 



a 



b 

~L0 



— + — + ao - 



(72) 
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The coefficients are functions of 5i, 82, A, and the inte- 
gration constant Cx necessary to match the z — ^ cx) and 
the z ^ series expansions. One can choose various 
combinations of those parameters that make some of the 
coefficients in the series expansion vanish. In particular, 
the choice 



62 
Si 



167r3/2 



16CA7r3/2' 
leads to the largest simplification 

79 77 (37r2 - 207rA + lOA^) 



(73) 
(74) 



Kret(z) 



4z4 12(7r-2A)Az2 
-2437r4 - SGOOtt^A + 233007r2A2 



A{z 



96(7r-2A)2A2 
-89607rA3 + 30740A4 
^ 96(^ - 2A)2A2 

The function A{z) near the origin behaves as 
TTZ^ (^-V^il + CxSi) + VXS2) 



(75) 



A3/2,52 



(76) 



_ z*{tt - 2A)(-48^3/2(^ ^ p^^^) ^ ^ 327162)) 

32A5/2J2 

With the conditions on 62 and 5i discussed above, this 
function becomes 



A{z) 



7r(7r-2A)2z6 



3\ ^8 



4A3 

TT (IStt^ - 827r2A + 1307rA2 - 52A' 
24A4 

The function B{z) near the origin behaves as 
52(7r-2A)2z4 



(77) 



Biz) 



+ CxSi) + VX52) A3/2 
a\z'^ 



(78) 



32 



tt(^-^{1 + CxSi) + VX52) A5/2 



where a\ is an involved function of A, C, Si and 62- 
Again, with the conditions on 62 and 5i discussed above, 
this function becomes 



B{z) 



(67r2 + 20A(A - 27r)) 
37rA - 6A2 



0{z') (79) 



3. Bom-Infeld in \ = n /2 

In this case, the strongest divergence of the 
Kretschmann goes as ~ z~20 _|_ 0{z~^^) + . . .. If we im- 
pose the following condition 



S2=V2{l + C^6i) , 



(80) 



then the leading order divergences are suppressed to ^ 
z~^'^ . Demanding that the coefficient of that term also 
vanishes, we find the condition 

327r 



<5i = 



V2 - 3271(7^! 



which leads to 



Kret(z) 



739 



64z4 15^2 



797 290787 



3200 



(81) 



(82) 



This result indicates that the modified gravitational ac- 
tion cannot soften the strength of the divergence below 
the limits of the standard BI model in GR. 
The function A(z) near the origin behaves as 



A{z) 



2(V2(l + (7^Ji)-52)^' 



(83) 



(5i - 16 



^(3^/2(1 + (7^(51) -2^2 



327r52 

With the conditions on 62 and 8\ discussed above, this 
function becomes 



A{z) 



9z 



10 



8 



+ 0(z"). 



(84) 



The function B{z) near the origin behaves as 
9<52^* 



B{z) 



V2{l + C^5i)-S2 



+ 0{z''). (85) 



With the conditions on 62 and 61 discussed above, this 
function becomes 



B{z) « -16- 



512z2 7084z4 598016^6 



15 



225 



23625 



(86) 



D. Inner horizons 



Given the units that we are using, for astrophysical 
black holes with low charge to mass ratio (Si <C 1), the 
existence and location of the outer horizon will be only 
slightly modified as compared with the GR case. For this 
reason, in this work we focus on the analysis of the inner 
horizons. The effects of the Planck-scale modified dy- 
namics on microscopic black holes, whose outer horizon 
may be significantly affected by the internal structure 
near the singularity, will be explored elsewhere. 

Horizons are determined by the cuts between the func- 
1/2 

tions 1 + 6iG{z) and 52zfj^ . In Fig|4]we can see how 
these curves change as we modify A and Z^. For the GR 
case (A = 0), for instance, we see that changing Ip may 
lead to the existence or absence of one inner horizon. To 
have an inner horizon we need a small enough value of Z^. 
The structure of A = is smoothly deformed as A is in- 
creased. As depicted in FigHl the case A = 7r/2 also has 
an horizon for small Ifs, though it disappears for larger 
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l+i5iG[z] 



1.0 




0.5 


^^5^5^=2- 10" V, 








-0.5 




2 3 


4 


5 


-1.0 










-1.5 










-2.0 











FIG. 4. Representation of the curves 1 + lt\G{^z) (various 
colors) and 52zf^^ (blue) for the GR case (red), the \ — 7v/2 
case (green) and the A = vr case (orange) for the relations 
r, = 10" and If, = IQ-V,, IQ-^r,, 2 x 10"^,, 5 x 10"^,. 



l+5iG[z] 



0.6 



0.4 



0.2 




FIG. 5. Representation of the curve 1 + SiG{z) for — 
1.1 X Ip'^ , Ip = If'* , and la = 0.9 x Ip* (lower, central, 
and upper curves respectively for each color) for A = (blue), 
A = 7r/2 (green), and A = 1/10 (red). Note that /^"'(O) / 

;^"*(V2)/ir*(i/io). 



values. For values of A larger than tt/2, the function G{z) 
develops a minimum and a maximum near the origin, as 
we saw in Figs[T]and[21 which may lead to the formation 
of up to two inner horizons for some combinations of pairs 
(A, Ip). In this section we describe the structure of inner 
horizons for the various combinations of such pairs. 



1. Critical value for Ip 

Let us define If^'^ as the value of above which there is 
no inner horizon. A natural way to determine the value of 
Z^*"'* is to look at the GR solution. In this case, the inner 

horizon disappears when the function AI{z) = l + 5iG{z) 
vanishes at the origin. The value of Ip that fulfills this 
condition turns out to satisfy the relation 1 -\-5iCgr — 0, 
being Cqr — G(0) (the integration constant needed to 
match with the z 3> 1 solution). This constraint gives 
^cnt jj-^ terms of Vq and rs as follows 



l+(5iG[z 



icrit 9 

la — r . 



(87) 



Remarkably, this relation not only holds for the GR case 
with C = Cgr but it does also hold for C — Cjl and the 
general case C = C\ with A < 7r/2 (see Fig (5]). 

The structure corresponding to A < 7r/2 is standard 
since only one inner horizon may arise if any. For A > 
7r/2 one can also find up to two inner horizons. The 
critical length in this case sets the value of Ip that 
first shows two inner horizons, being the innermost one 
located at z = (see FigIS]). For larger values oi Ip, 
the innermost horizon moves away from the origin and 
approaches the second inner horizon. Eventually, the two 
inner horizons coincide and for larger values oilp no inner 
horizon arises. This behaviour with A is explained by 
the fact that for A < 7r/2 the slope of the function 1 + 
5iG{z) near the center is positive, while when A > 7r/2 it 




FIG. 6. Representation of the curve 1 + 5iG{z) for Ip = \.lx 
Ip* (red), Ip = Ip* (green), and Ip = 0.9 x (orange) for 
A = 37r/2. Note that the minima occur at the same location 
Zmin — 0.7414, where dG/dz = 0. 



becomes negative, thus allowing up to two cuts with the 

1/2 

curve S2zfJ when I 



^ jcrit 



In this sense, the special 
case X — tt/2 represents the transition value between 
these regions and corresponding to a vanishing slope of 
the curve 1 + 5iG{z) at the center, while A = (the GR 
case) becomes an upper limit for this slope. Note that 
the value of the function 1 + SiG{0) does not depend on 
A, but just on the value of Ip as compared to the one of 

icrit 
Lp . 



2. Critical value for Si 

The previous discussion can be reinterpreted in a dif- 
ferent way. First we should note that our definition of 
G{z) is missing its dependence on the parameter A and, 
therefore, we should use the alternative notation G{z; A). 
This function G{z; A) can be expanded near the origin as 
G{z; X) = C\ + axz^ + bxz'^ + Therefore, the func- 
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FIG. 7. Evolution of M{zl^i„; Si) as Si is changed for A — 
7r/2, 37r/2, Stt (green, orange, and blue respectively). The x- 



axis represents the variable Si, which goes from 5i = to 5i 
in each case. Note that for A = 71/2 we have Si = S^'* . 



tion M{z; 5i) = l + SiG{z; A) at the origin takes the value 

AI (0; 6i) = l+SiC\. For those values of A for which C\ < 

0, the special condition M(0;^i) = picks out a special 

value of Si , vifhich we may denote as = -l/Cx. With 

this notation, we have G(z; A) = — tx + a.\z^ + bxz'^ + . . . 

°i 

, which explains why those configurations with Si > Si 
have M(0;(5i) = 1 - Si/S^ < and those with Si < S^ 
have M{0;6i) > 0. This notation also makes clear why 
the plots presented in Fig|S] for curves with different A 
but the same Si hit the vertical axis at the same point. 



3. Characterization of double inner horizons 

We have just seen that if M(0; 6i) < then we must 
have an inner horizon. However, M(0) > does not 
necessarily mean that we do not have an inner horizon. 
In this situation we may have none (if A < 7r/2) or two 
(for some values of <5i if A > 7r/2) or even one degener- 
ate (extreme). Now we will detail the conditions for the 
existence of two inner horizons for a given A. 

Consider FiglS] The minimum of those curves are lo- 
cated at the same position z^m- The functional depen- 
dence of z^j„ with A can be obtained analytically but the 
resulting expression is too long to be written here and not 
too much illuminating anyhow. Nonetheless, what really 
matters is that we know exactly where is located 

for a given A. The next step is to study the evolution of 
M(z4j„; ^i) = H-(5iG(z4mi ^) ^1 is changed (see Fig. 



ED. 



When the curve Af(z,^jj„; (5i) cuts the horizontal axis. 



-^i^min'^i *) = then that value of Si'* represents a 
degenerate double horizon. In other words, for Si — Si 
we have two inner horizons, being the innermost located 
at z = 0, where M{0;S^) = 0; for S^'* < Si < S^ we 
have two inner horizons; for Si = S^'* the two inner hori- 
zons converge at the same point; and for Si < S^'* there 



are no inner horizons. The curves appearing in FiglT] 
are segments of straight lines with an end located at 
{5^, M{z}^^„;S^),S^) and the other at (0, 1). This means 
that those lines can be expressed as follows 



1-M(zi,„;^^) 
6t 



Si . 



From this expression it is clear that Si'* 
0, which leads to Si'* — — 



y 



is defined by 
-. Note that 



l-M{z^.^;S^ 

both Si and M(z;^„;(Sf) can be computed numerically 
for each A. 



V. STUDY OF f{R) = R- R^/Rp 



In this case the equation (j35]) can be manipulated to 
allow for a more transparent interpretation and analysis. 
Taking into account the following identity 



1 



„3/2 
JR 



Jr + -j^fR.. 



1 



dz 



we can express ([55]) as follows 



4 V 47ry Iz/^/^ 



(89) 



(90) 



' R 



This expression can also be written as 



d 
dz 



M 



47r 



1 



1/2 dz 



(91) 

The advantage of this expression is not clear yet. How- 
ever, for the model f{R) = R — R^/Rp in which fa van- 
ishes at some small z, we will see that the right-hand side 
diverges as 1 / y^z — z+ (where z+ is defined below) and, 
therefore, its integral is finite. This means that the di- 
vergence of M(z) is totally described by the second term 
on the left-hand side. We thus have fully under control 
the divergent piece of the mass function. 

We now determine the series expansion of the right- 
hand side of (|9ip near the divergence. The divergence 
occurs when Jr — 0, i.e., at a curvature R ~ Rp/2. For 
this reason, we write this function as follows (A = A/vr) 



(92) 
(93) 
(94) 





-((T+ 

a 


— (T 


a = 




1 


<y± = 


h 


1 


1 


Vi + 


4A 


^± 


2A2 





1 + A~X±{1 + 2~X) 



(95) 
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FIG. 8. Representation of M{z) for GR (blue), and \ = 3n/2 
for f{R) ^ R + R'/Rp (green) and f{R) = R - R^/Rp 
(red). The curves represent the cases 5i = 5^^'^'^ ^"^ {1 ± 2 x 
10"^). Note that GR and f{R) = R - R^/Rp only present 
one horizon, whereas f{R) = R + R? /Rp may develop up to 
two horizons. 




FIG. 9. Representation of M{z) for GR (blue), and A = 37r/2 
for f{R) ^ R + R^/Rp (green) and f{R) = R - R^/Rp 
(red). The curves represent the cases Si — ^i^" and 5i — 
5i^"{l - 3 X 10"^). In these cases the GR solution 5i = 5^^° 



has the inner horizon at z 
horizon. The case f{R) = 
horizon as long as 5i > 0. 



= and Si < 5i has no inner 
R — R? / Rp develops an inner 



With this notation, it is easy to see that fn vanishes at 
z = For A ^ 1 it is also easy to see that z+ ~ 
A^/'^, where ttA = k'^/3^/Rp — Ip/lp and, therefore, z+ ^ 
Ip/lp. 

Note that (t_ — ^^^^^ that as A — > oo we find 

CT+ = 1+ whereas cr_ = 1^. This means that (t+ > 1 
always. It will also be useful to express A in terms of z+ 
as follows 



4(1 + 4) 



1 + + 2^4(1 + 4) 



(96) 



of ((9T|) near z = z_(. + e as follows 



RHS' 



47r 1 + zj 



(97) 



A/2 



16%/27r (1 + ^z+ + 2zl + 2zl^l + z* 



where e = — a. This shows that the integral of this 
function is finite, as claimed before. The function M{z) is 
plotted in Figsl8]and[9] As we see there, an inner horizon 
always exists, regardless of the model parameters when 
the constraint > is satisfied. 

From this analysis we conclude that for the model 
f{R) = R— / Rp the mass function diverges as M{z) ^ 
1 / •y/z"^^^+ as we approach z"*". 



The internal re- 
gion always has an inner horizon as long as > (which 
is the only physically meaningful situation). 

Given the expressions and (|97l) . the structure of 
the solution near the singular point z = z+ is completely 
known up to an integration constant coming from the 
integration of ([W]) . 



A. Metric components 

Using the same decomposition of the mass function 
M = 1 + SiG{z;X) as in the f{R) = R + R^ / Rp case, 
near the critical point z+ the metric functions A{z) and 
B{z) can be approximated as follows 



Ol 

<r2 



a2 

£3/2 



as 



04 



age 



1/2 



A{Z) = ^ + ^TTTT H h + flR 

= - + 4?2 + ^3 + b.e'l' + 

The form of the coefficients can be obtained recursively, 
being the leading terms given by 



.(98) 
(99) 



ai 



(-l + g+)5/V(5i 

167r (1 + 0-+) 3/2,52 
IGtt (-1 + a+){lj_a+y% 
-l + alSi 



(100) 
(101) 



B. Kretschmann scalar 

To study the degree of divergence of the Kretschmann 
scalar, we will consider the leading order contribution 
near z+, i.e., when e — 0, of the functions A{z) and 
B{z). For clarity, it is useful to replace the infinitesimal 
e — cr+ — CT by its corresponding expression in terms of 
C = z — z-(- . The correspondence can be easily found and 
is given by 



2C 



2{a\ - 1)5/4 



Using this expression, we can expand the right-hand side 



(z - z+) 



(102) 
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We thus have the following leading-order expressions as 
z — >■ z+ 

^(.)«^^ + + ... (103) 

S(z)«— ^ + ... (104) 
Inserting these expressions into (j59p . we find 

Kret{z) w ^- — ^ + . . . 

(105) 

This result is extremely interesting. It means that even 
though the metric components are strongly divergent as 
z+ is approached, the curvature at that point diverges 
just as ~ 1/(2: — z+)2, which is milder than in the case 
of GR and f{R) = R + R^/Rp where it goes as l/z^. 
This sharp difference in the intensity of the divergence is 
entirely due to the shift of the singularity from z = to 
z ~ z+, because in this case terms such as (—1 + / z'^ 
are no longer divergent. The divergent piece is just due 
to the last term of ((5^ . 



VI. CONCLUSIONS 

We have studied electrically charged black holes within 
the context of Palatini f{R) — R ± R^ / Rp gravities, by 
considering Born-Infeld electrodynamics. This choice of 
nonlinear electrodynamics as the matter source was cho- 
sen in order to break the tracelessness condition of the 
energy-momentum matter tensor, thus allowing to ob- 
tain solutions different from those of GR. In this work we 
have focused, in particular, in black holes whose charge 
to mass ratio is small (r^/rs ^ 1). These black holes 
only deviate from their GR counterparts near the center, 



and lead to several relevant modifications: i) new hori- 
zons may arise: in the f{R) — R + R^ /Rp case there 
may be up to two inner horizons (which have been fully 
characterized), having also the possibility of being de- 
generate (extreme), ii) in the f{R) — R — R^/Rp case 
the corresponding black holes show a similar structure in 
terms of horizons to that of the Reissner- Nordstrom solu- 
tion of Einstein-Maxwell equations in GR, but the radial 
coordinate cannot be extended below a minimum value 
r+, iii) the strength of the central curvature singularity 
may be softened due to Palatini /(i?) gravity: from r^^ 
of the Reissner-Nordstrom solution, and the of BI in 
GR (or r^* if the combination of the model and black 
hole parameters are properly chosen), we go to the 
in Palatini f{R) = R + R^ / Rp gravity but to (r - r+)-2 
in Palatini f {R) = i? - R^ / Rp gravity. 

Our results show that Palatini f{R) theories are able 
to severely modify the internal structure of electrically 
charged black holes. In the f{R) — R~ R^/Rp case, the 
softening of the divergence by the shift of the central sin- 
gularity to a finite point motivates the study of extensions 
of Palatini f(R) theories and also other NED models to 
better understand how the interplay between modified 
gravity and modified matter sources at very high curva- 
tures may affect this internal black hole structure and 
ameliorate its singularities, an issue to be explored in 
future works. 
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